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Abstract. Numerical simulations of hot accretion flows have shown that the mass accretion 
rate decreases with decreasing radius. Two models have been proposed to explain this result. 
In the adiabatic inflow-outflow solution (ADIOS), it is thought to be due to the loss of gas in 
outflows. In the convection-dominated accretion flow (CDAF) model, it is explained as because 
that the gas is locked in convective eddies. In this paper we use hydrodynamical (HD) and 
magnetohydrodynamical (MHD) simulations to investigate which one is physical. We calculate 
and compare various properties of inflow (gas with an inward velocity) and outflow (gas with an 
outward velocity). Systematic and significant differences are found. For example, for HD flows, 
the temperature of outflow is higher than inflow; while for MHD flows, the specific angular 
momentum of outflow is much higher than inflow. We have also analyzed the convective stability 
of MHD accretion flow and found that they are stable. These results suggest that systematic 
inward and outward motion must exist, i.e., the ADIOS model is favored. The different properties 
of inflow and outflow also suggest that the mechanisms of producing outflow in HD and MHD 
flows are buoyancy associated with the convection and the centrifugal force associated with the 
angular momentum transport mediated by the magnetic field, respectively. The latter mechanism 
is similar to the Blandford & Payne mechanism but no large-scale open magnetic field is required 
here. Possible observational applications are briefly discussed. 
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1. Introduction 

In the early analytical studies of hot accretion ^ow such as ADAFs, it was assumed 
that the mass accretion rate is a constant of radius. Correspondingly, the density follows 
p(r) oc r^'^l'^ . Later many numerical simulations have been performed, including both 
hydrodynamical (HD) and magneto- hydrodynamical (MHD) ones (e.g., Igumenshchev & 
Abramowicz 1999; Stone, Pringle & Begelman 1999; Stone & Pringle 2001). One of the 
most important findings is that the inflow and outflow rates (defined as the mass fiux of 
fiow with an inward and outward radial velocity) decrease with decreasing radius. These 
simulations, no matter they are HD and MHD ones, indicate that the profiles of inflow 
rate are quite similar, which can be described by, 

Afi„(r) = 27rr2 / pmin(w„ 0) sin^df? = Mi„(rout) f — ) , (l-l) 
h \routJ 

with s ~ 0.5 — 1 (see a review in Yuan, Wu & Bu 2012). Correspondingly, the radial 
profile of density becomes fiatter, p{r) oc r^P with p ~ 1 — 0.5. The consistency between 
HD and MHD results is surprising given that the mechanisms of producing the inflow 
rate proflle in HD and MHD cases are different, as we will illustrate in §2.3. But such a 
result is what predicted in the recent work of Begelman (2012). We want to emphasize 
that these theoretical results have been conflrmed by observations to, e.g., Sgr A* and 
NGC 3115 (see Yuan, Wu k Bu 2012 for details). 

A natural question is then: what is the nature of the inward decrease of the inflow rate? 
Two models have been proposed. One is the adiabatic inflow-outflow solution (ADIOS; 
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Blandford & Begelman 1999; 2004; Begelman 2012). In this model the inward decrease of 
inflow rate is because of the mass loss in the outflow launched at every radius. The second 
model is the convection-dominated accretion flows (CDAFs; Narayan, Igumcnshchev, & 
Abramowicz 2000; Quataert & Gruzinov 2000). This model is based on the assumption 
that the hot accretion flows, both HD and MHD, are convectively unstable. In this 
scenario, the inward decrease of mass accretion rate is because that with accretion, more 
and more fluid is locked in convective eddies operating circular motion. 

Answering this question is the first main of the work by Yuan, Bu & Wu (2012). In 
that work, we have investigated this problem from the following two aspects using HD 
and MHD numerical simulations. We first compare the various properties of both inflow 
and outflow, including the radial and rotational velocities, temperature, and Bernoulli 
parameter. If the CDAF scenario is correct, i.e, the motion of the flow is dominated 
by convective turbulence, and correspondingly inflow and outflow rates are simply due 
to turbulent fluctuation, we should expect that the properties of inflow and outflow are 
almost the same. As we will see, however, we find that the properties of infiow and outflow 
are systematically and significantly different (§2.1). The second aspect is to directly 
analyze the convective stability of an MHD accretion flow. The HD hot accretion flow is 
convectively unstable. However, there has been a debate on the convective stability of an 
MHD flow. Yuan, Bu & Wu (2012) used the simulation data to analyze the convective 
stability of MHD accretion flows and finds that the flow is stable (§2.2). Based on the 
above results, they have concluded that the decrease of accretion rate for both HD and 
MHD fiows is not because of convection, but systematic outflow. The origin of outflow 
and the observational applications are also discussed in §2.3 and §3, respectively. 

2. HD and MHD Numerical Simulations: ADIOS or CDAF? 

We have pcrfornwxi both HD (Model A) and MHD (Model B) simulations of two- 
dimensional axisymmetric accretion flows around black holes rising the ZEUS code (Yuan, 
Bu & Wu 2012). The initial condition of the two models is identical, which is a rotating 
torus with constant angular momentum. The difference is that magnetic field is included 
in Model B, which is confined to the interior of the torus. 

2.1. Different properties of inflow and outflow 

We have calculated the angle- integrated and mass flux-weighted value of some quantities. 
The results arc shown in Figs. 1-3. From the figures, we can clearly see that the properties 
of inflow and outflow are significantly different. For example, we see from Fig. 1 that for 
Model A, the temperature of outflow is significantly higher than that of inflow. From Fig. 
2 we see that for both Model A and B, the radial profile of the radial velocity (in unit of 
Keplerian velocity) is quite different for inflow and outflow, with the former increasing 
rapidly inward while the latter keeping roughly constant. From Fig. 3 we see that for 
Model B, the outflow has much larger angular momentum than the inflow. The difference 
between Model A and B also suggests that the mechanisms of producing outflow in the 
HD and MHD cases are different, as we will discuss in §2.3. This indicates that the inflow 
and outflow arc not simply due to convective (for Model A) and MHD (for Model B) 
turbulence. Rather, they are systematic inflowing and outflowing motion. 

2.2. The convective stability of MHD accretion flows 

We have analyzed the convective stability of an MHD accretion flow based on the H0iland 
criteria. We found that the accretion flow is convectively stable. This result indicates that 
CDAF model at least can't be applied to MHD flows. 
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Figure 1. The radial distribution of flux- weighted temperature for Model A (left) and B 
(right). The solid and dotted lines are for inflow and outflow, respectively. 




Figure 2. The radial distribution of the flux-weighted radial velocity for Model A (left) and B 
(right). The solid and dotted lines are for inflow and outflow, respectively. 




Figure 3. The radial distribution of the specific angular momentum / = v^rsinO for Model 
A (left) and B (right). The solid and dotted lines are for inflow and outflow, respectively. The 
dashed line denotes the Keplerian angular momentum at the equatorial plane. 



2.3. The Outflow Scenario and the Origin of Outflow in HD and MHD Flows 

The overall scenario suggested by the results mentioned above is that the inward de- 
crease of the inflow and outflow rates is because of the mass lost in outflow. Outflows 
occur throughout the accretion flow from almost any radius and they must escape out 
of the outer boundary of the accretion flow. But note that the existence of systematic 
outflow does not exclude the existence of turbulence. A snapshot figure of velocity vec- 
tors which indicate the coexistence of both systematic inward and outward motion and 
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turbulence is presented in Yuan, Bu & Wu (2012). A similar conclusion is reached in Li, 
Ostriker & Sunayev (2012), but Narayan et al. (2012) present a different point of view. 

Origin of outflow in HD flows was originally suggested to be due to the positive 
Bernoulli parameter of ADAFs. However, many numerical simulations, such as Stone, 
Pringle & Begelman (1999), found the existence of outflow, while the Bernoulli parameter 
is negative. On the other hand, from Fig. 1, we see that the temperature of inflow is 
significantly lower than that of outflow. This result, combined with the fact that HD 
accretion flows are convectively unstable, suggest that outflow is produced by buoyancy. 

Origin of outflow in MHD flows must not be because of the buoyancy since ac- 
cretion flow is convectively stable (§2.2). We see from Fig. 3 that the specific angular 
momentum of outflow is very close to the Keplerian value, much larger than that of inflow. 
This suggests that the outflow is produced by the centrifugal force. Consider two fluid 
elements located at two different radii in a differential rotating accretion flow. Magnetic 
stress transports the angular momentum from the inner fluid element to the outer one. 
Once the angular momentum of the outer element reaches nearly Keplerian value, the 
centrifugal force will be able to make the fluid element turn around and throw it outward. 
This mechanism is different from the Blandford & Payne mechanism since here we don't 
need a large-scale open magnetic field. In fact, we find that the magnetic field in both 
the accretion fiow and coronal region is tangled. We therefore call it a "micro-Blandford 
& Payne" mechanism. 

3. Possible Applications 

Many important problems remain to be probed. It is even not very clear in what 
condition the outfiow can reach infinity, although they are able to escape out of the 
outer boundary of the accretion fiow; and what is the exact values of outflow rate and 
terminal velocity. But still it is useful to discuss some possible observational applications 
which in turn can constrain theoretical models and supply some valuable clues. One 
is the Fermi bubble detected in the Galactic center. Our calculations indicate that the 
bubble can be readily inflated by the outfiow from the ADAF (Mou et al. in preparation) . 
Another example is the winds widely detected in AGNs and black hole X-ray binaries 
with various luminosities. It has been shown that in at least some sources a magnetic 
mechanism is required. People usually invoke Blandford & Payne mechanism. However, 
it is still unclear how to realize the required large-scale open magnetic field in accretion 
fiow. The micro-Blandford & Payne mechanism is then a promising alternative. The 
readers are referred to Yuan, Bu & Wu (2012) for details. 
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